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ABSTRACT 

A comparison of X-ray observations of the Cassiopeia A supernova remnant taken in 2000, 2002, 
and 2004 with the Chandra ACIS-S3 reveals the presence of several small scale features (< 10") which 
exhibit significant intensity changes over a 4 year time frame. Here we report on the variability of 
six features, four of which show count rate increases from ~ 10% to over 90%, and two which show 
decreases of ^ 30% - 40%. While extracted 1-4.5 keV X-ray spectra do not reveal gross changes 
in emission line strengths, spectral fits using non-equilibrium ionization, metal-rich plasma models 
indicate increased or decreased electron temperatures for features showing increasing or decreasing 
count rates, respectively. Based on the observed count rate changes and the assumption that the 
freely expanding ejecta has a velocity of ~ 5000 km s~^ at the reverse shock front, we estimate the 
unshocked ejecta to have spatial scale variations of 0.02 — 0.03 pc, which is consistent with the X-ray 
emitting ejecta belonging to a more diffuse component of the supernova ejecta than that seen in the 
optically emitting ejecta, which have spatial scales ~ 10~^ pc. 

Subject headings: ISM: individual (Cassiopeia A) - supernova remnants - ISM: dynamics 
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1. INTRODUCTION 

Cassiopeia A (Cas A) is currently the youngest 
known Galactic supernova remnant (SNR) with an es- 
timated explo sion date no earlier than around 1670 
AD ( Thorstensen et al.ll2003: iFesen et a l. 2006a). The 
remnant consists of optical, infrared, and X-ray emis- 
sion arranged in a shell roughly 4' diameter (~ 4 
pc at 3.41q '^ kpc; iReed et all 11995ft consisting of 
undiluted supernova (SN) ejecta rich in O, Si, Ar, 
Ca, and Fe heate d by the re mnant's 2± 30 km 
s~ reverse shock llChevalier fc Kirshn e^ Il978[ IT979t 
Douvion et al.'1999';'Huehes et al. 2000; Willineale et al.' 
2Q03 ; Hwane & Lamine 2003; Lamine & Hwane 2003; 
Morse et al.H2004|) . 



Viewed in X-rays, the remnant's brightest emission is 
concentrated in a 30" thick, ~ 110" radius emission ring 
dominated by thermal plasma arising from shock heated 
SN debris. Farther out beyond this emission ring, lies 
fainter filamentary X-ray emission (radius ~ 150") as- 
sociated with the current position of the remnant's « 
6000 km s~^ forward blast wave as it moves through the 
surrounding circumstellar material ijGotthelf et alJl2001t 
IPeLancv & Rudnick 2003). 

If the remnant's expanding SN ejecta were uniform in 
density and chemical composition. X-ray emission from 
the shocked ejecta would appear as a uniform ring of 
emission with little spectral variation. What is actually 
observed, however, is a patchy, irregular ring of clumpy 
emission reflecting the inhomogeneous nature of the ex- 
panding SN debris as it interacts with the reverse shock 
front. Due both to strong ejecta clumping caused by ra- 
diative cooling instabilities and the turbulent mixing of 
different chemical layers of the progenitor star during the 
SN explosion, the resulting X-ray emission morphology is 
structurally and chemically complex ( Hughes et al. 200Qi: 

^ Smithsonian Astrophysical Observatory, 60 Garden St, Cam- 
bridge, MA 02138 

^ 6127 Wilder Laboratory, Physics & Astronomy Department, 
Dartmouth College, Hanover, NH 03755 



'Hwa ng et a'l]l2000|) . 

Although Cas A's reverse shock front has been located 
to arcsecond accuracy in a few regions through high- 
resolu tion optical images ijFesen et al.ll200H iMorse et alJ 
l2004j) , it has only been broadly localized through a sharp 
increase in X-ray emission in azimuthal averages of se- 
lected portions of the remnant (Gotthclf et al. 2001). On 
small spatial scales, evidence for the reverse shock inter- 
acting with the ejecta would be manifested in X-rays as 
a rise in emission as the SN ejecta is shocked and decel- 
erated. Post-shock cooling of the metal-rich ejecta could 
lead to equally dramatic localized emission declines. 

Here we present an analysis of archival Chandra images 
of Cas A which show a number of substantial, small- 
scale brightness changes between 2000 and 2004 which 
we interpret as marking regions of reverse shock passage. 
We report on spectral analyses for six specific regions 
and find some plasma temperature changes in the X-ray 
emitting ejecta. In §2 we present the relevant Chandra 
observations and subsequent data reduction. In §3 we 
discuss the results from our analysis, and summarize our 
conclusions in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

Cas A was observed with the ACIS-S3 chip as part of 
the Chandra Guest Observer program for 50 ksec both 
on 30 January 2000 and 6 February 2002, and for 1 
Msec as a Chan dra Very-Large Project during the first 
half of 2004 fseelHwang et alil200nH Gotthelf et al.1l200ll 
iDeLanev fc Rudnickll2003l fo r details r egarding the 2000 
and 2002 observations, and IHwang et al.i i20M for the 
2004 observations.). The 0"492 CCD pixel scale of ACIS 
undersamples the telescope's ~ 0"5 resolution. We re- 
processed these data using the latest version (3.2) of the 
Chandra Calibration Database (CalDB). This reprocess- 
ing focused on recalculating the aspect solution for better 
positional accuracy. Version 3.2 of the CalDB was also 
used in the generation of detector response matrices used 
in the spectral fitting, which is discussed below. 

Both direct visual comparisons of the three epoch im- 
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ages and image differences and subtractions were used 
to identify several small scale regions where the detected 
flux appeared to vary with time. Six selected regions 
showing significant flux variations in the total ACIS band 
(0.3 - 10.0 keV) are marked in Figure ^ with enlarge- 
ments shown for the years 2000, 2002, and 2004 in Fig- 
ure [3 The coordinates, sizes of the extracted flux re- 
gions, and count rates are given in Tabled The quoted 
rates for each epoch are normalized to the observed 
count rate of the remnant's central X- ray point source 
jTananbaum. .1999.; .Aschenbach 19991 iPavlov fc ZavlinI 

Spectra were extracted for the six regions listed in Ta- 
ble n Background subtraction was done by selecting 
regions close to the knots. The background spectra were 
binned such that each detector channel contained a min- 
imum of 10 counts. The background spectra were then 
fit to models appropriate for a thermal plasma in non- 
equilibrium ionization. More specifically, since we are 
investigating regions in the interior of the SNR, the back- 
ground contains contributions from swept up circumstel- 
lar material as well as diffuse emission from shocked 
supernova ejecta around the regions of interest. This 
method of background subtraction is preferred to using 
ACIS Blank-Sky images since scattered X-rays from the 
SNR are the dominant source of background here. Fur- 
thermore, while the source spectra were found to vary, 
the background spectra comes from the diffuse emission 
around the knots and does not vary. 

In order to minimize the effects of Chandra's az- 
imuthally varying point spread function w hich are man- 
ifest in the combined 2004 observations ijHwang et al.l 
120^, we limited ourselves to using a 50 ksec subset 
(ObsID 5196; 8 February 2004) of the 1 Msec observation 
when extracting the 2004 spectra, since the observational 
parameters for these data closely match those of the 2000 
(ObsID 114) and 2002 (ObsID 1952) observations. The 
extracted source and background spectra were fit using 
version 11.2 of the X-ray spectral fitting package XSPEC. 

The six regions of interest are relatively small in area 
(< 30 arcsec^) with total observed counts of 5000 - 20,000 
for the 0.3 ~ 10 keV energy a range for each of the three 
epoch observations. We grouped the extracted spectra 
such that there was a minimum of 10 counts in any de- 
tector channel. Moreover, since the goal of this analysis 
was not concerned with looking for specific abundance 
variations in the ejecta . we chose to limit ou r spectral 
fitting to NEI models l|Borkowski et al.ll200l1) contain- 
ing mainly Si, S, and Ar (and Mg and Ca, when there 
were prominent emissions visible at 1.2 and 3.8 keV, 
respectively) and limited the fits to between 1.0 and 4.5 
keV. For all models, the ratio (Si/Si©) was fixed at 1, 
while the relative abundances of the other included ele- 
ments were allowed to vary. Finally, we chose to freeze 
the abundances of all other modeled elements (H, He, 
C, N, O, Ne, Ca, Fe, and Ni) at zero. We justify this 
by noting that these elements do not have emission lines 
in the energy band we are interested in. We do note, 
however, that by not including these elements, we may 
under-predict the contribution to the fiux from the un- 
derlying continuum. In Table |21 we list the fitted plasma 
temperature for each epoch and each region, and show 
the data and the resulting fits in Figures |3HS| 

All regions, with the exception of Region 5, are well 



modeled as a thermal spectrum. Region 5 appears to 
possess purely nonthermal emission and is likely filamen- 
tary emission arising from the remnant's forward shock 
front viewed in tangent lik e several other inter ior pro- 
jected filamentary features ('Hu ghes et al. ''2000: Region 
'D') and common in the r emnant's southwestern quad- 
rant l|DeLanev et al.ll2004(l . 

3. RESULTS AND DISCUSSION 

3.1. Count Rate Variations 

A comparison of 2000 - 2004 Chandra images of Cas 
A shows fairly steady X-ray emission flux from the rem- 
nant's shock heated ejecta on both large and small spa- 
tial scales. However, there are a few small features in 
the main emission ring that do show significant changes 
in brightness over this time interval. Most of these show 
increases rather than decreases although this may be due 
in part to a detection bias against finding changes in dif- 
fuse patches of emission, especially if they are faint. 

As can be seen qualitatively in Figure Inland quantita- 
tively in Table HI between 2000.1 and 2004.1 four smaU 
regions in Cas A's X-ray emission structure exhibited sig- 
nificant increases in intensity while two regions showed a 
decrease. Specifically, Regions 1, 3, and 4 each show in- 
creases of « 10%-20% between 2000 and 2004, while Re- 
gion 2 went from near invisibility in 2000 to being among 
the brightest features in the remnant, with a near dou- 
bhng of its count rate between 2002 and 2004. In 2000, 
the count rate in this small region is due mainly to the 
diffuse emission in the area around and toward the knot 
in Region 2. As can seen in the 2000 and 2002 images, 
a new and relatively bright knot emerges here, becoming 
twice as bright just two years later. 

We note that fine-scale flux changes in the remnant 
are not necessarily monotonic. Region 3, for example, 
showed both an increase and a decrease over four years; 
nearly a 40% increase from 2000 to 2002 but then a 14% 
decrease from 2002 to 2004. The short time span of 
data taken in early 2004, taken as part of the IMsec 
Very Large Program on Cas A, does not permit us to 
investigate whether this decrease continued after Febru- 
ary 2004, but future observations could investigate this 
question. 

Most of these flux changes occur on relatively small 
spatial scales, of order a few arcseconds. For example, 
the white box in Figure 2 centered on the emerging bright 
spot in Region 2 is just 3" x 3" and yet encloses most of 
the observed count rate change. Similarly, the majority 
of the brightness change for the northeastern limb emis- 
sion knot shown in Region 1 occurs along a segment of 
the knot's structure extending 1" - 3" to the northwest. 

In contrast to the brightening seen for Regions 1-4, 
the X-ray emission detected in Regions 5 and 6 show de- 
creases in brightness. Shown in the bottom two rows of 
Figure |21 the count rates in these regions drop signifi- 
cantly from 2002 to 2004; a ~ 40% decrease for Region 
5 and ^ 20% in Region 6. For Region 5, the measured 
40% change in flux represents an average over several 
small, individual features some of which exhibited large 
and small changes in brightness. 

While it has been known for some time that a layer 
of contamination is forming o yer the A CIS detector thus 
reducing its overall efficiency ijPlucinskv et al. .2003) , the 
expected flux reduction from this effect is not enough to 
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explain the overall decreases we see here. Nonetheless, to 
avoid uncertainties due to detector sensitivity degrada- 
tion, we list in Table 2 measured count rates normalized 
to the observed count rate for the remnant's central X- 
ray point source. Conversion of the observed count rates 
into intensity changes can be estimated using the X-ray 
point source ob served Fx = 8 x 10~^'^ erg c m~^ (for 

0. 6 — 6 keV; see lFesen. Pavlov. &: Sanwaill2006l and refs. 
therein). 

3.2. Spectral Variations 

We performed spectral fits for each region to investi- 
gate whether rapid changes in the observed count rate 
also reflected changes in spectral properties of the emis- 
sion features. In general, we find that the modeled spec- 
tra are relatively constant from year to year (cf. , Region 

1, Fig. 01 left, and Regions 4, Fig.^ right), but there are 
exceptions. For example. Region 2 (Fig. O right) which 
showed the most extreme flux brightening, shows some 
evidence for the emergence of helium-like silicon and sul- 
fur emission starting in the 2002 spectrum. From our 
spectral flts for this region, we also found that kT in- 
creased in value from ~ 1.0 - 1.5 keV between 2002 and 
2004, combined with an increase in n^t from 4.6 - 53 x 
10^*^ cm~^ s over the same time period. Such large in- 
creases in flux and ionization timescale suggest that the 
ejecta in this location has indeed recently encountered 
t he remnant's rev e rse sh ock front. 

lAnderson et al.l l|1994|) predict that ejecta knots which 
are strongly decelerated by the reverse shock will show 
both a marked increase in their thermal emission (like 
that seen here), coupled with an increase in synchrotron 
emission. We have compared 6 cm VLA radio observa- 
tions (~ 0'.'4 resolution) taken in April 2000 and again 
in April 2001 and note a qualitative increase in the radio 
emission around Region 2 ( Dclanev 2006) which is com- 
parable to the increase in X-ray emission. This increase 
in the radio and X-ray flux i s con sistent with the inter- 
pretation of lAnderson et alJ l)1994|) . However, we do not 
find any corresponding changes in the radio emission for 
the other five regions we present here, contrary to what is 
e xpected if we interpret t hese flux changes in the context 
of lAnderson et alJ 1)19941) . In particular. Region 1, which 
like Region 2, also shows a change in kT (from 0.8 - ~ 
1.1 keV) between 2000 and 2004 , does not exhibit a cor- 
responding increase in the 6 cm radio emission ( Dclanev 

nam . 

On the other hand, model fits for Region 4, a large dif- 
fuse patch along the eastern rim and which showed about 
the same level of X-ray flux changes as seen for Region 
1, exhibited nearly a constant electron temperature over 
the four year period. Moreover, Region 3 (Fig. 0] left) 
did not exhibit the same sort of strong, monotonic kT 
changes as the other three regions. However, based on 
the large change in the modeled kT between 2000 and 
2002 (Table|2Il from ~ 0.8-1.5 keV, along with the emer- 
gence of a broadening near helium- like sulfur (tentatively 
identifled here as Si Ly/3) suggests that the ejecta here 
may be undergoing the same sort of rapid ionization seen 
in Region 2. Another possible spectral change seen for 
Region 3 is the apparent broadening of the silicon lines in 
the 2002 and 2004 spectra. Combined with the increase 
in kT, this suggests that the ejecta is in the process of 
ionizing up to the hydrogen-like state for silicon. 



From its spectrum. Region 4 (Fig. H right) appears 
to be nearly in ionization equilibrium, and as noted 
above (Table |2J) the fitted temperature stays nearly con- 
stant between 2000 and 2004. Modeling suggests that its 
spectrum can be equally well modeled by either a NEI 
plasma with an ionization age ~ 10^'^ cm~^ s, or by a 
collision al ionization model, such as a Raymond-Smith 
plasma ijRavmond fc Smithlll977f) . The high fitted ion- 
ization age is puzzling because it implies electron densi- 
ties of ~ 10^ cm""^; densities which are as high as those 
seen in the optical knots. However, visual inspection of 
the spectrum for Region 4 (Fig. 01 right) shows that it is 
clearly different from regions which do not exhibit ion- 
iz ation equilib r ium ( i.e.. Regions 1 and 3 for instance). 

orse et ahl |2004) proposed that Cas A's ejecta con- 
sists of cool, optically bright emission knots embedded in 
a warm, tenuous ejecta medium which is responsible for 
the observed X-ray emission. This general description is 
supported by our model fits (Table EJ. Were the ejecta 
completely uniform, then a single temperature, single 
ionization age would adequately describe the spectra in 
all regions at each epochs. However, since we see that 
the fitted temperature rises in Regions 1-3, the rate at 
which the ejecta is being shocked in these regions is not 
constant due to density inhomogeneities in the preshock 
material. 

While the four regions which showed increases in in- 
tensity appear to be associated with ejecta which has 
recently crossed the reverse shock. Regions 5 and 6 show 
flux decreases and therefore do not indicate reverse shock 
locations. The interior projected thin filament seen in 
Region 5 has a non-thermal spectrum and is likely as- 
sociat ed with the forward and not the reverse shock 
front l|DeLanev et all 120041) . The observed chan ges m 
its brightness probably reflects changes in path length 
along interconnecting tangents of the wave front of the 
forward shock as it interacts with the local circumstellar 
medium on the near- facing side of the remnant. 

In contrast, emission seen in Region 6 (Fig. O right) 
which also showed a decrease in brightness, may repre- 
sent ejecta cooling in the post-shock flow following re- 
verse shock passage. Its projected location in the rem- 
nant is consistent with this identification. While it lies at 
approximately the same radial distance from the center 
of expansion as the brightening Regions 1-4, it lies just 
behind (eastward) a patch of optically bright emission 
ejecta along the remnant's southwestern limb. If this 
feature is indeed fading due to cooling, this would im- 
ply it was shocked much earlier than in the other regions 
in the area and would be consistent with its location in- 
side the steady X-ray and optical emission features found 
farther out from the remnant's center. 

3.3. Scale Length of X-ray Emitting Knots 

Rapid intensity changes in small regions throughout 
the remnant offer some insight into the structure of the 
lower density component of Cas A's SN ejecta. Abrupt 
increases (or decreases) in the observed count rate 
strongly suggest that the X-ray emitting ejecta is non- 
uniform on fairly small scales. This conclusion has long 
been clear from optical observations (Fcscn et al. 2001| 
but not from X-ray observations which are sensitive to 
a much lower density, higher temperature component of 
the ejecta but which likely dominate the remnant's over- 
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all mass budget Jyink et al.lll99()HWillingale eiUMM 
iLaming fc Hwangll200.lD . The velocity of the freely ex- 
panding ejecta clumps at the specific regions we investi- 
gated is not known, but a rough estimate can be made 
based on the distance of the reverse shock from the cen- 
ter of expansion. The reverse shock is located Vrs ~ 95" 
(~ 1.6 pc at a distance of 3.4 kpc) from the center of 
expansion l)Gotthelf et alJl2001|) and unshocked, freely 
expanding ejecta at this distance have been expanding 
for tage — 330 yr. Therefore, we estimate that the un- 
shocked ejecta velocity is no more than v^j ~ frs/tage ~ 
5000 km s^^ 

The observed intensity for the small features we stud- 
ied show significant changes on time scales of At ~ 2-4 
yrs. This suggests a scale length for spatial variations 
in the ejecta no larger than the distance covered by a 
~ 3000 km s~^ reverse shock moving through ~ 5000 
km s~^ expanding ejecta over the course of a few years; 
i.e.. At X Vej ~ 0.02 — 0.03 pc or about 1-2 arcsecond at 
3.4 kpc. This scale length is consistent with the actual 
sizes of the variable flux features we see in the remnant. 
In contrast, the fine-scale structure of optical knots i n 
Gas A are of order ~ 0'.'2 - 0'.'3 in size l|Fesen et al.ll200lD . 
roughly an order of magnitude smaller than the spatial 
scale for variations seen in the clumpy component of the 
X-ray bright ejecta. 

4. CONCLUSIONS 

With the exception of the rapidly evolving remnant as- 
socia ted with SN 1987A l)Biirrows et a l. 2000':' Park et all 
l2004j) , no remnant has been previously reported to show 
significant flux changes across portions of its X-ray emit- 
ting structure. This is due in part to the limited number 
of high-resolution images available on young SNRs cover- 
ing a significant time span, and the relatively small num- 
ber of very young Galactic or LMG/SMG SNRs known. 

Inspection of X-ray data taken of the Gassiopeia A su- 
pernova remnant with Chandra AGIS-S3 in 2000, 2002, 
and 2004 reveals the presence of several small scale fea- 
tures (< 10") which exhibit significant intensity changes 
over this 4 year time frame. Here we described six such 
features four of which had count rate increases from 
10% to over 90%, while two others showed decreases of 
~ 30% - 40%. Whereas the majority of the six variable 
flux features exhibited no gross spectral changes in the 
1-4 keV energy band during the four year period, one 
region which had the greatest increase in brightness also 
showed hints of slight emission line changes as well as an 
increase in its 6 cm radio emission. Spectral fits using 
non-equilibrium ionization, metal-rich plasma models in- 



dicate increased or decreased electron temperatures for 
features showing increasing or decreasing count rates, re- 
spectively. 

A fading emission filament projected near the rem- 
nant's interior (Region 5) showed a completely nonther- 
mal spectrum and thus may be associated with the for- 
ward shock front. It decreasing intensity may be the re- 
sult of changing path lengths along the edges of a wavy 
forward shock as it interacts with the local, inhomo- 
geneous circumstellar medium. Based on the observed 
count rate changes and the assumption of freely expand- 
ing 5000 km s~^ ejecta at the reverse shock front, we 
estimate the unshocked ejecta to have spatial scale vari- 
ations of 0.02 — 0.03 pc, consistent with the X-ray emit- 
ting ejecta belonging to a more diffuse component of the 
supernova ejecta. 

In general, bright X-ray and optical knots are not co- 
incident. This is not unexpected given the differences in 
their temperatures and densities. Moreover, the sizes of 
the X-ray knots 1" — 5") appear to be an order of 
magnitude larger than the optical knots 0"2 — 0"3; 
iFesen et al.ll200l|) . However, the remnant's optical knots 
can also appear and brighten on timescales as short as a 
few years, not unlike the timescales we observed for these 
X-ray knots. Given the newly observed short term vari- 
ability seen here in the X-ray, future X-ray and optical 
observations of selected knots and filaments might pro- 
vide valuable information regarding the dynamical and 
radiative evolution of supernova ejecta knots over a range 
of densities, sizes, and postshock temperatures. 

Gassiopeia A is currently the youngest and one of the 
most studied Galactic SNRs. The small yet noticeable 
changes we found combined with future X-ray observa- 
tions of this bright X-ray remnant may help shed light 
not only on the remnant's near-term X-ray evolution, but 
also help us better understand the fine-scale structure, 
dynamics and chemistry of debris from core-collapse su- 
pernovae in general. 
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TABLE 1 
ACIS-S3 Count Rates 



Region 


RA (J2000) 


Dec (J2000) 


Extraction 


Normalized Count Rate^ 


Count Rate Cliange 




(h:m:s.s) 


(O , 


Region 


(2000.1) 


(2002.1) 


(2004.1'' ) 


(A2004-2000) 


1 


23:23:34.0 


+58:49:56.2 


5" X 4" 


1.84+0.02 


2.04+0.01 


2.09+0.01 


+14% 


2 


23:23:24.4 


+58:47:13.2 


3" X 3" 


0.48+0.02 


0.57+0.01 


0.93+0.01 


+93% 


3 


23:23:13.4 


+58:48:10.0 


5" X 3" 


0.93+0.01 


1.28+0.02 


1.10+0.02 


+18% 


4 


23:23:11.2 


+58:48:53.2 


5" X 5" 


2.06+0.02 


2.24+0.02 


2.51+0.01 


+22% 


5 


23:23:26.9 


+58:48:19.0 


3" X 10" 


1.56+0.02 


1.21+0.02 


0.95+0.02 


-40% 


6 


23:23:18.0 


+58:48:37.0 


4" X 4" 


1.04+0.01 


0.83+0.01 


0.71+0.02 


-31% 



^Values listed are detected counts s ^ normalized to Cas A's X-ray point source rate of 0.12+0.01 cts s ^. 

^Epoch 2004 count rates are calculated from Chandra Observation ID 5196, which has the same roll angle as the 2000.1 and 2002.1 observations. 
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TABLE 2 

Ti':.Mi'i' u.vn 



HE \AHI.\TIOXS 



Region 


2000.1 


kT (keV) 
2002.1 


2004.1 


1 

2 
3 
4 
6 


94+"-»2 
'^•''^-0.03 


n Sfi+0 09 

i.uy_o 28 
1 43+0.07 
-■-■^■^-0.10 

Q2+0-04 

n f;o-|-0.02 


, .c+0.17 

1 45+0.08 
^•*^-0.12 
94+0 02 

n p:c+0.02 
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Fig. 1. — A Chandra image of Cas A in 2004 (Observation ID 5196). The marked regions correspond to tliose listed in Table [Tl and 
shown in Figure 1^ 
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Fig. 2. — Close-up of the six regions listed in Table [Tl for each epoch. The boxes correspond to the spectral extraction regions used in 
the fits in Table 1^ and shown in Figures 
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Fig. 4.— Left: Spectral fits to Region 3 for 2000, 2002 and 2004. Right: Same as left, for Region 4. 
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Fig. 5.— Left: Spectral fits to Region 5 for 2000, 2002, and 2004. Right: Same as left for Region 6. 



